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ABSTRACT 



We show here that the rotation period data in open clusters allow the empirical determination 
of an expression for the rate of loss of angular momentum from cool stars on the main sequence. 
One significant component of the expression, the dependence on rotation rate, persists from prior 
work; others do not. The expression has a bifurcation, as before, that corresponds to an observed 
bifurcation in the rotation periods of coeval open cluster stars. The dual dependencies of this 
loss rate on stellar mass are captured by two functions, f(B — V) and T(B — V), that can be 
determined from the rotation period observations. Equivalent masses and other [U BVRIJHK] 
colors are provided in Table 1. Dimensional considerations, and a comparison with appropriate 
calculated quantities suggest interpretations for / and T, both of which appear to be related 
closely (but differently) to the calculated convective turnover timescale, r c , in cool stars. This 
identification enables us to write down symmetrical expressions for the angular momentum loss 
rate and the deceleration of cool stars, and also to revive the convective turnover timescale as a 
vital connection between stellar rotation and stellar activity physics. 

Subject headings: convection - open clusters and associations: general - stars: activity - stars: 
evolution - stars: late-type - stars: rotation 



1. Introduction 

Cool stars are known to lose angular momentum and spin down over time. However, a detailed under- 
standing of this loss has yet to be achieved. Measurements from the ground and space, coupled with theory, 
give us a good idea of the rate of angular momentum loss for the present-day Sun, but extensions to stars 
of different masses or those of other ages, particularly very young stars, are problematical in various ways. 
This paper proposes a path toward making such extensions. 

In the same paper that proposed the existence of the solar wind, Parker (1958) noted that this wind 
would cause a 'retardation of solar rotation' on a timescale similar to its lifetime. Weber & Davis (1967) 
elaborated on this particular effect and showed that the associated angular momentum loss rate is equivalent 
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to assuming corotation of the solar wind out from the surface to the Alfvcnic radius, ta- Schatzman 
(1962) had previously suggested that the solar case could be generalized to other stars, linked the associated 
angular momentum loss to the presence of a surface convection zone, and provided a formula for the loss 
rate in terms of a star's angular velocity and certain other quantities. However, it was not clear how to 
calculate or measure these quantities. A similar criticism could be leveled against the expression for angular 
momentum loss provided by Mestel (1968). The associated viewpoint that the angular momentum loss rate, 
dJ/dt oc Qr\M where f2 and M are the star's angular velocity and mass loss rate respectively, ultimately 
runs into the difficulty of providing credible calculations or measurements of both ta and M for the Sun and 
cool stars of other masses. Alternatively, writing d,] /dt oc BqR 4 £1/va, where Bo, R, and va are the surface 
magnetic field, stellar radius, and wind velocity at the Alfven surface, requires deriving or measuring Bq and 
va as a function of a star's mass and other properties. See Collier Cameron & Li (1994) for a discussion of 
the issues that arise when one takes this viewpoint. 

Another way to approach this problem is empirical. Kraft (1967) measured a decline of stellar rota- 
tion velocities with age, strengthening the solar-stellar analogy, and Schatzman's work provided a natural 
framework for interpreting the break in the Kraft (1967) curve of rotation velocities against stellar mass. 
The work of Skumanich (1972) provided more direct information about the angular momentum loss rates 
for solar-mass stars when he noted that their measured rotation velocities, wsini, decline with age, t, as 

vsini oc t~ ' 5 (1) 



where fsini represents an average over coeval cluster stars. This observation implies that the rate of loss of 
angular momentum for (solar-mass) stars of constant structure (i.e. moment of inertia) obeys the relationship 

According to this viewpoint, if the observations can provide the acceleration, dQ/dt as a function of all 
relevant variables, then dJ/dt immediately follows if the associated moment of inertia can be inferred. 

A hybrid approach helps overcome certain weaknesses of both viewpoints. Endal & Sofia (1981) used the 
wind-loss formulation of Belcher & MacGregor (1976), in turn based on the work of Weber & Davis (1967), 
to understand the rotational history of the Sun, but the resulting rotational evolution (their Figure 4), 
including certain assumptions about internal transport of angular momentum, did not closely match the 
empirical results of Skumanich (1972). At that time, the security of the Skumanich result was not assured, 
and its applicability was restricted to stars of near-solar mass. Indeed, Mestel (1984) investigated other loss 
relationships than dJ/dt oc — £1 3 and made significant steps toward elucidating what they implied. Kawaler 
(1988) extended Mestel's work in certain significant ways, identifying dependencies in addition to in the 
expression for angular momentum loss, and proposed an expression that could be implemented easily in 
models of rotating cool stars of varying mass. The relationship, combining theoretical and observational 
considerations, is 

^ = -K w n 1+ ^ 3 ( A)2-» ( M 14 )l-2n/3 ( (3) 

at riQ Mq 

where R, M, and M are the stellar radius, mass, and mass loss rate respectively, n parameterizes the geometry 
of the magnetic field, a is the exponent that relates the surface magnetic field, Bo to £1 via Bq oc £l a , and 
K w is a calibration constant chosen to ensure that a solar model attains the solar rotation rate at solar age. 
n can vary from 3/7 for a dipolar field to 2 for a purely radial field. If one assumes a linear dynamo, then 
a=l, and n must then be set equal to 1.5 to reproduce the observed Skumanich spin-down, vsini oc i~ 0,5 . 
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This choice of n kills the mass loss term, of course. Consequently, the above expression usually simplifies to 

ft =-^<^»°' 5 (^)"° s - w 

eliminating the mass loss rate, M. This expression, a hybrid of the two viewpoints discussed above, has 
been routinely used to drain angular momentum from rotational stellar models constructed using YREC, 
the Yale Rotating Stellar Evolution Code (e.g., Pinsonneault et al. 1989), and also is used in other stellar 
models (e.g., Bouvier et al. 1997). A didactic account of the above developments, set in a broader context 
than here, may be found in Chapters 13 and 21 of Maeder (2009). 

We note that some studies of the rotational evolution of stars do not explicitly provide an expression for 
angular momentum loss (e.g., MacGregor & Brenner 1991; Armitage & Clarke 1996). This is usual in cases 
where the magnetic field and associated quantities (and sometimes their evolution) are themselves modeled 
numerically. The mass loss rate is usually an input parameter. A recent example is the work of Matt & 
Pudritz (2008), which provides details of such an approach. Some models include the possible effects of disks 
on the pre-main-sequence (e.g., Collier Cameron & Campbell 1993, Keppens et al. 1995, Sills et al. 2000, 
Barnes et al. 2001). But here we are concerned only with the main-sequence where any possible effect of 
pre-main-sequence disks has abated. We also note that cataclysmic variable (CV) research uses an expression 
derived by Rappaport et al. (1983) for the extra-gravitational loss of angular momentum. This expression 
has some similarity to the Kawaler (1988) expression, as noted by Andronov et al. (2003). However, CVs 
are beyond the scope of this paper. 

In the observational domain, rotational data for stars in open clusters had revealed (and further ob- 
servations confirmed) the presence of the so-called ultra-fast rotators (UFRs) in the Pleiades (van Leeuwen 
& Alphenaar 1982). These observations implied that the angular momentum loss rate for such stars was 
much lower than that calculated using Kawaler's expression. Chaboyer et al. (1995), following MacGregor 
& Brenner (1991), therefore suggested modifying this relationship at high rotation rates to facilitate the 
theoretical modeling of such stars. This idea became known as "saturation" (e.g., Stauffer 1994) and was 
roughly parallel to a similar phenomenon in soft X-rays, (although the exact connection between the two 
was not clear then). Chaboyer et al. (1995) suggested the loss relationship 

*L -K (—\0-5(^_\-o.5 j ^ forO<fi crit , 

dt Kw[ R Q ] [ m g > \ni it n forO>tt crit , [) 

where S! crit is a constant. (A way of achieving a similar effect in the context of Weber-Davis type wind models 
is presented in Collier Cameron & Li 1994.) Furthermore, Barnes & Sofia (1996) showed that the UFRs 
could not be modeled by the original Kawaler expression, with its f2 3 dependence, regardless of any possible 
pre-main-sequence spin-up, the bloated state of T Tauri stars being considered a possible reservoir of angular 
momentum suited to explain the origin of the UFRs. As a result, subsequent work, e.g., Krishnamurthi et 
al. (1997) and Barnes et al. (2001) routinely used the newer Chaboyer et al. (1995) loss prescription. It was 
hoped that fi C rit could simply be set equal to some constant threshold angular velocity for all relevant cool 
stars. 

However, the accumulation of additional data suggested that a single constant value of f2 cr i t for all cool 
stars was inadequate (Barnes & Sofia 1996; Krishnamurthi et al. 1997). Consequently, Krishnamurthi et al. 
(1997) advocated a scaling 

^crit = ^crit©— (6) 
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where t c is the convective turnover timescale. (Prior to this, and in the context of using a Weber-Davis 
type wind to understand the rotational evolution of stars, Collier Cameron & Li (1994) had also suggested 
a scaling involving the convective turnover timescale, in this case setting the star's surface magnetic field, 
B according to B oc T c fl.) 

The Ohio State University group (Sills et al. 2000; Andronov et al. 2003; Denissenkov ct al. 2010) 
absorb the solar radius and mass from the Kawaler/Chaboyer formula into the wind constant, K w , and 
write: 

dJ „ [R [ n 3 for ft < Q crit , 



(7) 

dt 1M[!1^0 forO>ft crit , 

but ft cr it is now variable as described immediately above, or varied piece-wise to match open cluster obser- 
vations. For example, Sills et al. 2000 used the Krishnamurthi et al. (1997) scaling for ft cr ;t for masses 
down to O.6M0, and set ft C rit manually for each modeled stellar mass below this. Stars of very low mass are 
beyond the scope of this work, but we note that such a tuning effectively modifies f2 C rit from a constant into 
a new and arbitrarily modifiable function. Making f2 cr i t mass dependent implies, of course, that the y/R/M 
term in the expression for dJ/dt is simply not capturing the mass dependence of the angular momentum 
loss, as it was originally intended to. Clearly, there is a problem. 

On the other hand, the bifurcation in the angular momentum loss, as proposed by MacGregor & Brenner 
(1991), encapsulated in the relationship of Chaboyer et al. (1995), and as subsequently used in modeling 
is important in retrospect, because it does capture an essential feature of the observations. Indeed, since 
then, a steadily growing rotation period database has allowed the identification of distinct fast (C-) and 
slow (I) sequences in color-period diagrams of open cluster stars, as initially proposed by Barnes (2003). 
This C /I classification has been confirmed by extensive rotation period observations in M 35 (Meibom et al. 
2009), M37 (Hartman et al. 2009) and M50 (Irwin ct al. 2009), the first including a decade-long radial 
velocity survey for cluster membership and multiplicity. Scholz & Eisloffel (2007) include a discussion of this 
bifurcation vis-a-vis rotation periods in Praesepe/Hyades, and most recently, data in Hartman et al. (2010) 
clearly display this bifurcation in a large and uniform rotation period study of the Pleiades. 

This observed bifurcation ties in well with the bifurcation in the angular momentum loss expression 
proposed by Chaboyer et al. (1995), although not to the mass dependence used there, which is the same for 
both kinds of stars. We shall see below that the mass dependence is actually different for the two sequences, 
so that it is impossible for the same expression to describe both. 

Barnes (2003) also suggested interpretations of the observed shapes of the C- and I sequences, based on 
theoretical considerations, and a unifying scenario (hereafter called the Cgl scenario) proposed in that work 
for the rotational evolution of cool stars. In particular, he suggested that stars initially are fast rotators on 
the C sequence, where the inner radiative and outer convection zones are largely decoupled, so that the mass 
dependence of this sequence is specified by (the reciprocal of) the moment of inertia of the outer convection 
zone. The observed transition of stars from the C- to the I sequence was proposed to be coincident with a 
change in the mass dependence, which was itself proposed to change from that of the outer convection zone 
alone to that of the entire star, and also to be coincident with the onset of an interface dynamo. Thus the 
mass dependence for the I sequence was suggested to be dependent on (the reciprocal of the square root of) 
the moment of inertia of the whole star. 

This work begins by considering and checking whether these proposed dependencies work. We show that 
they do not (Section 2). We then show that the observations themselves might be queried to provide dJ/dt 
directly in terms of observed quantities (Section 3). Section 4 proposes an interpretation of the relevant 
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observed quantities in terms of the convective turnover timescale. The relationship with stellar activity 
physics is pointed out in Section 5, and the conclusions are stated immediately thereafter. (The next paper 
in this series combines the C- and I-type behaviors identified in Section 4 into a simple nonlinear model for 
the rotational evolution of stars and explores the consequences for gyrochronology.) 

2. Inadequacy of the moment of inertia proposal 

We begin by negating a proposal made by Barnes (2003), that the mass dependence of the rotation 
periods in open clusters can be simply attributed to the moments of inertia of either the whole star or that 
of the surface convection zone. 



2.1. I sequence 

Kawaler (1989) was the first to note that beyond an age of a few hundred million years, as exemplified 
by the 600Myr-old Hyades open cluster, a deterministic relationship between rotation period, color and 
age could be derived, and inverted to provide a star's age. Younger clusters have a more complicated 
morphology, initially parsed into fast/C- and slow/I sequence stars by Barnes (2003). He proposed that the 
latter, I sequence stars are describable by Pj = f(B — V) x g(t) where f(B — V) and g(t) are empirically 
determinable functions of the B — V color and age, t, respectively. In that paper, the functional forms used 
were f(B — V) = \J{B — V — 0.5) — 0.15(B — V— 0.5) and g(t) = \ft. These functions have been subsequently 
re-determined, most notably by Meibom et al. (2009), based on a very large study of both rotation periods 
and membership in the open cluster M35. They determined that f(B - V) = 0.77(1? — V — 0.47) 55 , and 
for definiteness, we will use this latter form in this papeiQ. (Transformation to mass or other colors in the 
set [UBVRIJHK] can be accomplished using Table 1.) 

Barnes (2003) suggested identifying f(B — V) with l/y 7 ?*, where /* is the moment of inertia of the star. 
We have calculated /* using the latest version of YREC, and display it in Figure 1 as a function of mass and 
of B — V color for a series of 500 Myr models of solar composition. This age was selected to ensure that 
all lower-than-solar-mass stellar models of interest have passed through the pre-main-sequence phase and 
arrived on the main-sequence, while higher-than-solar-mass models have not evolved off the main-sequence. 
To the precision of this work, further main-sequence evolution does not have an appreciable effect. The 
numerical values and associated [UBVRIJHK] colors using both Green et al. (1987) and Lejeune et al. 
(1997), Lejeune et al. (1998) color transformations are provided in Table 1. 

We compare f(B — V) with \j\fh, in Figure 2. The curves are normalized to agree for the solar case. 
It is clear that the agreement is not good. One could argue that the normalization makes it hard to tell how 
poor the fit is for the cooler stars. The key part of the disagreement, however, is that does not drop 

off sufficiently fast for the warmer stars. This disagreement would only get worse if the normalization were 
changed. 



1 In a later comparison, we will also show the Mamajek & Hillenbrand (2008) determination of f(B — V). We also acknowledge 
that the index n in g(t) = t n is usually found to exceed 0.5 slightly (e.g., Barnes 2007, Collier Cameron et al. 2009, James et 
al. 2010). 
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Fig. 1. — Left: The moment of inertia of the entire star, /*, against stellar mass, in the range 1.5 — 0.15A/q . (Note 
the inverted mass axis.) Right: The same against B — V color. As expected, I* increases monotonically with stellar 
mass. (In this, and subsequent figures, we display the Lejeune et al. (1997, 1998) B — V color; other colors in the set 
[UBVRIJHK], including those with the Green et al. (1987) color transformations are provided in Table 1.) 

2.2. C sequence 

Barnes (2003) proposed that the C sequence of faster rotating stars can be described by Pc — Poe t / T< -- B ~ v \ 
where Pc denotes the rotation periods of the stars in question, t is the age, and T(B — V) is the relevant 
spin-down timescak@. The open cluster observations indicate that T(B — V) is short for stars bluer than the 
Sun, tending toward zero at mid-F. On the other side, redward of the Sun, T(B — V) is known to increase. 
The main requirement is a function that dives to zero at B — V = 0.47, the x-intercept of f(B — V), as 
determined by Meibom et al. (2009). The most reliable empirical determinations of T(B — V) for G and K 
stars are also by Meibom et al. (2009), who find that T(B — V) = 60Myr for G stars (specifically, those 
with 0.6 < B - V < 0.8) and T(B - V) = 140 Myr for K stars, (with 0.8 < B - V < 1.3). We estimate 
uncertainties on these values of 5 Myr and 13 Myr respectively, and use these values to represent those for 
mean B — V colors of 0.7 and 1.05. For redder stars, we note that about a third of the Hartman et al. 
(2009) sample of stars from the 550 Myr-old open cluster M37 are on the C sequence. This implies that 
T(B - V) = 500±65Myr for early M stars, represented by a point at B - V = 1.45. T(B - V) is plotted 
using these discrete values in Figure 3. We also display a cubic spline fit to these data, and acknowledge 
that T(B — V) becomes increasingly uncertain with decreasing stellar mast[f|. Similar results - 50 — 100 Myr 
for 1M & , 0.5 - 1 Gyr for 0.5M Q , and several Gyr for M < O.3M - were quoted by Scholz et al. (2009) and 
references therein. 

Barnes (2003) suggested that T(B — V) oc I cz , where I cz is the moment of inertia of the convection zone. 
We have calculated I cz using YREC, and display it in Figure 4 as a function of mass and of B — V color for 
a series of 500 Myr models of solar composition. I cz drops at both ends as expected; it decreases for stars 
with masses lower than 0.8 Mq, and also for the warmer stars with thinning surface convection zones. The 



2 The choice of an exponential function in t is partly motivated by theoretical considerations, and it is possible that other 
functions of t are also suitable, as discussed below. T(B — V) is our main concern here. 

3 The uncertainties quoted above are based on fractional numbers of stars, and do not include other effects such as uncer- 
tainties in the assumed cluster isochrone ages. 
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Fig. 2. — Comparison of the empirical function f(B — V) for the Isequence (dashed blue), with l/i/I* (filled circles). 
The latter was proposed by Barnes (2003) to interpret the empirical dependence. (The /* curve is normalized for 
solar mass.) We conclude that the two curves are not as similar as was proposed in Barnes (2003), negating the 
corresponding identification. In particular, 1/yJ* does not drop steeply enough to model the observed rotation 
period decline among the warmer stars. 



associated numerical values and other colors are provided in Table 1. 

We compare T(B — V) with I cz in Figure 5. The two variables are in reasonable agreement for M > M Q 
but diverge badly for M < M Q . Whereas T(B — V) is a steadily increasing function of B — V (or decreasing 
stellar mass), I cz becomes relatively flat for K dwarfs, and indeed drops for M dwarfs, as might be expected. 
Consequently, we find that I cz is a poor match to T(B — V), and withdraw the corresponding suggestion 
made in Barnes (2003). 



3. Empirical approach 

It transpires that the observations themselves can be queried to provide an empirical loss rate, through 
a logical extension of prior work. This is accomplished here, considering first the I type rotators, and then 
the C type ones. 



3.1. Isequence 

We return to the observational basis that the rotation periods, Pj, of stars on the Isequence are de- 
scribable as a product of two separable functions, /, and g, of the star's mass or B — V color and age, t, 
respectively: 

Pi = f(B -V)x g(t). (8) 

This was shown by construction in Barnes (2003) and Barnes (2007), where both f(B — V) and g(t) were 
determined empirically. 

g(t) = Vt (9) 
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Fig. 3. — The C sequence spin-down timescale, T(B — V), for cool stars (in units of Myr) plotted against B — V 
color. The green line is a cubic spline fit to the data. T(B — V) increases steadily with B — V color, somewhat faster 
than linear for the coolest stars. The data points are discussed in the text. 



was found there to be applicable to a wide range of cool stars, a result tracing back to the work of Skumanich 
(1972), where the more restrictive case of solar-mass stars was considered. As in the prior section, we will 
use the Meibom et al. (2006) determination 

f(B — V) = 0.77(B — V — 0.47) 055 (10) 

over other determinations, including those of Barnes (2003) and Barnes (2007), because of the size and 
quality of the data set considered there. 

Changing variable from P] to via Pj — 2tt/U, equation (8) becomes 

^ = f(B-V)xg(t), (11) 

and upon differentiating it with respect to time, we get 

dfl fl 2 Tr .dg 



Using g(t) — y/i from Equation (9), we get 



Re- using Equation (11) to write \ft in terms of / and J7, we get 

which specifies the deceleration of the star in purely observational quantities. Note that (the square of) 
f(B — V) completely specifies the mass dependence. f(B — V) being zero for mid-F and earlier-type stars, 
the deceleration is correspondingly zero for these stars. f(B — V) increases steadily through late-F, G, K, 
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Fig. 4. — Left: The moment of inertia of the outer convection zone, J cz , against stellar mass. (Note the inverted 
mass axis.) Right: The same against B — V color. J cz drops at both ends as expected; to the left as surface convection 
zones thin out and vanish, and to the right, with diminishing stellar mass. 



and early M stars, so that the deceleration of these stars is correspondingly greater. This accounts for the 
shape of the I sequence in open clusters. It is also worth emphasizing that the open cluster data directly 
provide only the deceleration above, not the angular momentum loss rate, dJ/dt, derived below. 

Deriving the rate of angular momentum loss requires an additional step. By the Chain Rule, 

dJ dfl dl 

— = I— + n— (15) 
dt dt dt y ' 

so that simple substitution from Equation (14) above gives 

Barnes (2003) suggested, in agreement with results from helioseismology, that the relevant moment of inertia 
/, in this case is /*, the moment of inertia of the entire star. To a very good approximation, this is constant 
on the main-sequence, killing the final term in Equation (16), so that the final empirical expression for the 
rate of loss of angular momentum, dJ/dt, from I sequence stars becomes 



This expression has the appealing feature that f(B — V) can be straight-forwardly determined from 
open cluster rotation period data, provided that the memory of the possibly complex initial conditions has 
been erased. We know this to be at least approximately true on the I sequence, because the Sun is on 
this sequence, and because helioseismic results indicate that the Sun is a solid-body rotator to first order. 
Consequently, expression (17) combines a relatively well-known and easily calculated function /», with the 
observationally well-determined function f(B — V) to state the mass dependence of angular momentum loss, 
independent of assumptions beyond those encapsulated in Equations (8) and (9) above. On the negative 
side, it has the undesirable feature of mixing calculated and observed quantities in one expression. It is also 
not explanatory, because the origin of f(B — V) is not yet specified. In a section below, we propose an 
identification for f{B — V). 
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Fig. 5. — The C sequence spin-down timescale T(B — V), in units of Myr (pluses; dashed green line), is compared 
with J cz , the moment of inertia of the surface convection zone in cool stars (unfilled circles; solid red line). J C z is 
normalized for solar mass. Although the agreement is good for M > Mq, partly because both I cz and T drop to 
zero at similar color, I cz and T(B — V) can be observed to diverge for lower mass stars, negating the corresponding 
identification proposed in Barnes (2003). 



We note that assuming / oc 1/y/TZ, as suggested in Barnes (2003), would simply replace i*/ 2 with 
a constant, and would not provide the correct mass dependence observed in the rotation period data, as 
shown above. Consequently, we here withdraw that suggestion. Also, expression (17) above can be compared 
directly with the f2 < fi C rit part of the prior loss rate expression (7) which is: 

*L = -K wX [^n 3 (is) 

dt V M y ' 

Comparing this with Equation (17) above shows that, apart from numerical constants, the difference in the 
loss relationships boils down to the difference between y/R/M and I*/ 2 . 



We plot both I*/ 2 and y/R/M against B — V color in Figure 6, so that a graphical comparison can 
be made. As before, the curves are normalized for a solar- mass model. Redward of B — V = 1.4, for very 
low mass stars, the physics is expected to be different, and neither relationship might be relevant. So, only 
the area blueward of this is relevant. Unsurprisingly, R/M, and consequently, y/R/M, is a relatively fiat 
function of B — V color for FGK stars, because the stellar radius tracks the mass. Consequently, the angular 
momentum loss rate for the prior prescription is essentially independent of spectral type, and any assumed 
initial distribution of rotation periods can be expected to retain this initial shape when subjected to this loss 
rate. Consequently, with the prior loss prescription, the observed shape of f{B — V) would have to be 
assumed as an initial condition. It would not arise as a natural consequence of the angular momentum loss 
prescription. 

In contrast, /*/ 2 is only somewhat flat for K stars (B — V ~ 1), where it attains its maximum (Figure6), 
and declines sharply blueward for G stars, dropping to zero blueward of late-F. This mass-dependence has 
been derived from the data, so it must be true within the uncertainties of the data themselves and the above 
assumptions. Consequently the loss rate is maximized at K stars, but declines steadily for G stars, dropping 
to zero at mid-F. Here, almost any non-pathological set of initial periods would progressively attain the shape 
of f(B — V) with the passage of time. 
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Fig. 6. — Comparison of the empirical dependence, I,f 2 (filled circles), proposed here for the Isequence, with the 
prior suggested color dependence of angular momentum loss, yf(R/M), (unfilled squares). The curves are normalized 
for the solar case. yjRjM is observed to be significantly different from the empirically determined function J*/ 2 . 



3.2. C sequence 

Here, we derive the second half of the angular momentum loss expression, for the C sequence, in a 
manner symmetrical to that for the Isequence above. We can tackle the C sequence in a manner similar to 
the Isequence, although without recourse to separability of the mass and age dependencies, by writing the 
rotation periods, Pc, of the C sequence stars as a function 



Pc = Pc(B-V,t) 



(19) 



of the B — V color and age, t, of the star. The open cluster observations show that the rotation periods of 
the C sequence stars have a roughly exponential behavior with agc@, suggesting that we write the periods of 
these stars in a form originally suggested by Barnes (2003): 

Pc = Poe t/TiB ~ v) (20) 

where Pa is the initial period for C sequence stars, and T(B — V) is the appropriate mass-dependent timescale 
for spin-down on the C sequence, as discussed in Section 2.2 earlier. The open cluster data currently provide 
T(B — V) only in discrete form. Figure 3 displays T(B — V) as a function of B — V color, along with a cubic 
spline fit. 



Switching variable from Pc to ft using Pc = 2tt/VL : equation (|20|) becomes 

HZI - p nP */T(B-v) 

n 



and on differentiating it with respect to time, we get 

dfi _ n 2 

~dt ~ 



t/T(B-V) 

2vr Q T(B-V)' 



(21) 



(22) 



4 An exponential is a compact way of expressing this dependence because it can be expanded in a power series in t/T and 
the coefficients determined. For example, for late-type stars with large T, the exponential will devolve into a function that is 
linear in t/T and will also lead to the same final result for dQ/dt. 
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Rcsubstituting from Equation (20) and again using Pc — 2tt/H, wc get 

dfi _ n 

~dt ~ ~T{B-V)' 



(23) 



which again specifies the deceleration of the star in purely observational quantities, this time for C sequence 
stars. Note that the mass dependence of the deceleration is specified completely by 1/T(B — V). Given 
the form of T(B — V) from Section 2.2, with T(B — V) increasing monotonically with B — V color, we 
see at once that as far as C sequence stars are concerned, the deceleration is greatest for late-F stars, and 
declines steadily as the mass decreases through G, K, and M stars. We note that the data provide only the 
deceleration above, not the angular momentum loss rate, dJ/dt, derived below. However, dJ/dt is required 
for comparison with prior work, and is potentially important in deriving conclusions about the mass loss 
rate and the Alfven Radius. 

Deriving the angular momentum loss rate requires additional steps. Applying the Chain Rule as before, 

dJ dSl dl c 



so that substitution from above gives 



dJ n t j_o dIc OK\ 

^ = -T{B—V) C+ ^f (5) 



We do not know yet what Ic, the relevant moment of inertia on the C sequence, is; but if we assume that it 
is not varying with time, we get as the final empirical expression for the rate of loss of angular momentum, 
dJ/dt, from C sequence stars 

As with the I sequence case, T(B — V) is derivable from the data, so that the rate of loss of angular 
momentum immediately follows provided that we can guess the appropriate moment of inertia Ic, for 
C sequence stars, and provided that this moment of inertia is not varying with time. (We will revisit this 
last assumption in a subsequent publication.) As with the I sequence, T(B — V) is empirical, and thus not 
explanatory, because its origin has not been specified. In the next section, we propose an identification for 
T(B - V). 

There is little choice about (a time-independent) Ic- It could either be the moment of inertia of the 
entire star, /*, or, if one followed the suggestion of Barnes (2003), the moment of inertia of only the outer 
convection zone, I cz . If Ic = I*, then dJ/dt would increase steeply with stellar mass, as shown in Figure 7 
(left panel, filled circles) , because the increase in 2* with stellar mass is amplified by the concomitant decrease 
in T(B — V). If Ic — Icz, then the mass dependence of dJ/dt would be far less steep with increasing stellar 
mass, as shown in Figure 7 (left panel, unfilled circles), and undefined blueward of B — V = 0.45, where both 
I cz and T(B — V) arc zero. (The two curves coincide for very low mass fully convective stars.) 

Expression (26) can be compared directly with the O > O cr i t part of the prior loss rate expression (7) 
which is: 

A t I D 

(27) 




We see that the difference between them boils down to the difference between I C /T(B — V) and ^/R/MQ 
The case of a constant O crit is displayed graphically in Figure 7 (left panel, squares), normalized to the 



crit' 



- 13 - 




Fig. 7. — Left: I*/T(B — V) (filled circles) and I C2 ,/T(B — V) (unfilled circles), both candidates for the empirical 
C sequence mass dependence, are compared with the prior \J R/M C sequence mass dependence (squares). While the 
first two decline steadily with increasing B — V color, y/R/M is flat. Right: The same functions 7*/T and I C z/T 
are compared with the Krishnamurthi et al. (1997) proposal that f2 C rit oc l/r c (instead of constant as at left). We 
see that the last is similar to I,/T(B — V) and I CZ /T(B — V), but more complicated. 



solar mass model with Ic — I*- We note that, as expected, it is a remarkably flat function of stellar 
mass, suggesting quite different expectations for the corresponding angular momentum loss rate than for the 
Ic = I C z and Ic = I* cases discussed above. Because of various issues discussed in their paper, Krishnamurthi 
et al. (1997) decided to abandon a constant Q C rit> and suggested a scaling where f2 crit oc l/r c . We display 
the result in Figure 7 (right panel, squares), normalized to the solar mass model with Ic — I*- We find that 
this scaling gives results similar to the Ic = I* or Ic — I C z cases above, which are both conceptually simpler. 



3.3. Combined expression 



To summarize the conclusions of this section, combining the two expressions above, the empirical accel- 
eration rate becomes 

dVl J n /r(s-y) for the C sequence 

~dt = ~ \ & f (B-V) for the T sequenC e ; ' ' >N ' 

where T(B — V) and f(B — V) are derivable from open cluster rotation period data. This leads immediately 
to the following expression for the angular momentum loss rate 

dJ _ _ U^{ Ic It{b-v) - dIc /dt} for the Csequence 

~dt ~ ~yi*i t f-(B-v)i^ 2 for the I sequence, ( ' 

and on further assuming that Ic , the relevant moment of inertia on the C sequence is constant in time, we 
get 

dJ \ nic /T(B-v) for the C sequence 

~dt ~ ~ 1 n3l *f 2 ( B ~V)/ s „ 2 for the Isequence. 
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4. Interpretation of T(B - V) and f(B-V) 

We now consider, in turn, possible interpretations for T(B — V) and f(B — V). We will suggest that 
both are related, albeit differently, to the convective turnover timescale in cool stars. 



4.1. T(B — V) 

We recall that T(B - V) is defined by 

Pc = P,e t ' T ^ v \ (31) 

where Pc represents the rotation periods of C sequence stars in open clusters, Pq is the appropriate initial 
period, and t is the age. T(B — V) is the appropriate mass-dependent timescale. Accordingly, T(B — V) 
must have dimensions of time. 

In Section 2.2, where it was first introduced, T(B — V) was discussed in detail. Here, we merely reiterate 
that the observations currently define it in terms of four data points. It drops to zero at B — V ^0.47, and 
increases steadily for lower mass stars, as shown in Figure 3. 

We are struck by the similarity between T(B — V) and the convective turnover timescale, r c , in cool 
stars, which we have calculated and list in Table 1 as a function of stellar mass and [U BVRIJHK] colors. 
r c and T(B— V) are displayed and compared in Figure 8 (left) and Figure 8 (right) on linear and logarithmic 
scales respectively. In this comparison, and the next, we use the global turnover timescale. However, the 
(somewhat shorter) local timescale at the base of the convection zone is also provided in Table 1. More 
details about the calculation can be found in Kim & Demarque (1996). 
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Fig. 8. — Left: r c (asterisks) and kcT(B-V) (pluses) are plotted against B—V color on a linear scale. Normalization 
at solar mass requires kc = 0.65d/Myr. A cubic spline fit to the T(B — V) data points (dashed green line) is also 
displayed. Right: The same comparison on a logarithmic scale, emphasizing the blue region where both functions 
decline sharply. We see that r c and kcT(B — V) are in reasonable agreement. 



The behaviors of r c and T(B — V) suggest the identification 



T(B-V) = £-, (32) 
kc 
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where r c is the convective turnover timescale, and kc is a (dimensionless) constant. We find that kc ~ 0.65 
d/Myr. Indeed, with this choice of kc, the two curves are in reasonable agreement. (At this early stage, it 
is not possible to comment on the dctcctability or importance of the slight non-monotonicity of t c for the 
lowest stellar masses.) 



4.2. f(B - V) 

The situation with f(B — V) is similar to that with T(B — V) above. Recall that f(B — V) is defined 
by the expression 

P l = f(B -V)x g(t), (33) 

where Pi represents the rotation periods of I sequence stars in open clusters, and f(B — V) and g{t) are 
separable functions, respectively, of the B — V color and age, t, of the stars. To a precision adequate for the 
present purposes, g(t) = \ft. Because Pj has the dimensions of time we infer that f(B — V) must have the 
dimensions of \/tirne, or that f 2 (B — V) must have the dimensions of time. 

The rotation periods of the I sequence stars decrease sharply for stars bluer than the Sun, approaching 
zero for B — V = 0.45. For stars redder than the Sun, the rotation periods increase all the way till at least 
0.5Mq. This behavior is again remarkably reminiscent of the behavior of the convective turnover timescale 
but dimensionally ought to be associated with its square root. We therefore propose the identification 

f(B-V) = ^, (34) 

where ki is a dimensionless constant, and r c is again the convective turnover timescale. The factor of 2 is 
added purely for later convenience. We find that fcj ~ 740 or 1340 Myr/d, where the smaller value is relevant 
to the f{B — V) determination of Meibom et al. (2009) while the larger value is relevant to the determination 
of Mamajek & Hillenbrand (2008). The functions 2r c and kif 2 (B — V) are compared graphically in Figure9 
and indeed they are in reasonable agreement. 



4.3. Symmetric combination 

The above interpretations of T(B — V) and of f(B — V) are appealingly symmetric. Although there 
are two timescales in the problem, captured by the dimensionless constants, kc and kj, the spin-down of 
both the C- and I sequences in open clusters appear to be connected by one underlying variable, r c , the 
convective turnover timescale in cool stars. As a result, it appears that the underlying variables in stellar 
rotation appear to be P, t, and r c , the last encapsulating the dependence on stellar mass. 

Consequently, placing the identifications above into Equation (28), we suggest the following expression 
for the deceleration of cool stars: 

dfi | n kc l T for the C sequence 

— = — < 3 (35) 
dt y Uir 2 k! for the I sequence 

or, on changing variable from O to P = 27r/f2, we get 

dP i kcP /r c for the C sequence 
\ T " /k,p for the I sequence 
which displays the symmetry of the problem explicitly. 
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B-V B-V 



Fig. 9. — Left: 2r c (asterisks) and kjf 2 (B — V) are plotted against B — V color on a linear scale. The dashed green 
curve is the Meibom et al. (2009) determination of f(B — V), while the dotted blue curve is that of Mamajek & 
Hillenbrand (2008). Normalization at solar mass requires ki = 740Myr/d and 1340Myr/d respectively. Right: The 
same comparison on a logarithmic scale, emphasizing the blue region where all the functions decline sharply. We see 
that 2r c and fe// 2 (B - V) are in reasonable agreement. 

5. Relationship with stellar activity physics 

In a classic paper, Noyes et al. (1984) proposed that the fractional chromospheric emission from cool 
stars was dependent, not just on the rotation period, P, but also on stellar mass. Inspired by theoretical work 
by Durney & Latour (1978), they identified the convective turnover timescale, t c (B — V), as the most likely 
variable to capture the dependence on stellar mass, and demonstrated that the fractional chromospheric 
emission of cool stars, R' HK , was more tightly dependent on the Rossby Number, Ro = P/t c , than on P 
alone. Soon thereafter, the transition-region ultraviolet emission and coronal X-Ray activity of cool stars 
were also shown to behave similarly (Vilhu 1984; Simon et al. 1985), and it is now generally accepted that 
stellar activity is associated more closely with the Rossby Number than with the rotation period alone. 

This work brings the rotational evolution of cool stars under the same umbrella by showing that dP/dt 
is a function of P/t c . Consequently, the rotation rate of a cool star appears to be a deterministic function 
of its age and its mass, the latter appearing through the convective turnover timescale, r c . 

6. Conclusions 

We have shown here that open cluster rotation period data allow the empirical determination of an 
expression for angular momentum loss from cool stars on the main-sequence. The relationship has a bi- 
furcation, as before, that corresponds to an observed bifurcation in the rotation periods of open cluster 
stars. 

One component of the relationship, the dependence on rotation rate, persists from prior work. The 
remaining dependence is captured by two observationally determinable functions of the mass or color of a star, 
which we write as f(B—V) and T(B — V) for the two tines of the relationship, respectively. Transformations 
to mass or other desired colors in the set [UBVRIJHK] can be accomplished using Table 1. 
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We derive an empirical expression for the angular momentum loss rate, which is 

dJ J n { Ic It — dIc I ' dt} for the C sequence 



and which simplifies to 



(37) 

dt in/,/ f or t ne I sequence 



d J \ nic It for the C sequence 

= — S 3 2 ("") 

dt \ n I 'i / S7T 2 for the I sequence. 



if the relevant moment of inertia, Ic , on the C sequence, is not time- varying. 

Both f(B — V) and T(B — V) appear to be related to the convective turnover timescale, t c in stars. We 
suggest the identifications 

T(B -V) = -j- (39) 

kc 

and 

/ 2 (B-^)-^, (40) 

where kc and ki are two dimensionless constants appropriate to the C- and I sequences, and are respectively 
approximately 0.65d/Myr and 740 or 1340Myr/d, depending on which of two particular forms one uses for 
f(B-V). 

Consequently our final expression for the angular momentum loss rate is 
dJ \ fJ { kcIc j Tc - dIc /dt} for the C sequence 



(41) 

dt y n I * Tc /iTT 2 ki for the I sequence 

and that for the deceleration of cool stars is 

dP ) cF I t c _ for the C sequence 
( Tc AiP f° r the I sequence 

which is both symmetric and dimensionally correct. We see that in this model, the evolution of the rotation 
period, P, depends only on the age, t, and the convective turnover timescale, r c , which encodes the mass 
dependence of rotation. 

Finally, we have pointed out that this model for the rotational evolution of stars makes a natural 
connection to stellar activity physics, where P/t c is the preferred independent variable. 

These ideas are developed further in following papers into a simple nonlinear model for the rotational 
evolution of cool stars from C- to I-type, and a companion paper will provide more general grids of related 
calculated quantities on and off the main-sequence. 
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Table 1. Calculated Properties of Solar Metallicity Main-sequence Stellar Models at 500 Myr, as used in the text 
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